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ABSTRACT 

Context. Messier 35 (NGC 2168) is an important young nearby cluster. Its age, richness and relative proximity make it a privileged 
target for stellar evolution studies. The Kepler K2 mission recently observed it and provided high accuracy photometric time series 
of a large number of sources in this area of the sky. Identifying the cluster’s members is therefore of high importance to optimize the 
interpretation and analysis of the Kepler K2 data. 

Aims. We aim at identifying the cluster's members by deriving membership probabilities for the sources within 1° of the cluster’s 
center, going further away than equivalent previous studies. 

Methods. We measure accurate proper motions and multi-wavelength (optical and near-infrared) photometry using ground based 
archival images of the cluster. We use these measurements to compute membership probabilities. The list of candidate members from 
|Barrado y Navascues et al.] ( |2001[ > is used as training set to identify the cluster’s locus in a multi-dimensional space made of proper 
motions, luminosities and colors. 

Results. The final catalog includes 338 892 sources with multi-wavelength photometry. Approximately half (194452) were detected 
at more than two epochs and we measured their proper motion and used it to derive membership probability. A total of 4 349 candidate 
members with membership probabilities greater than 50% are found in this sample in the luminosity range between 10 and 22 mag. 
The slow proper motion of the cluster and the overlap of its sequence with the field and background sequences in almost all color- 
magnitude and color-color diagrams complicate the analysis and the contamination level is expected to be significant. Our study 
nevertheless provides a coherent and quantitative membership analysis of Messier 35 based on a large fraction of the best ground- 
based data sets obtained over the past 18 years. As such, it represents a valuable input for follow-up studies using in particular the 
Kepler K2 photometric time series. 
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1. Introduction 


Different galactic environments, it is assumed, produce distinct 
populations during the star formation episodes. From very loose 
and young associations such as Taurus-Auriga, doomed to be 
dispersed very soon in the galactic gravitational well, to the very 
old halo population in massive globular clusters, there is a large 
gamut in ages, compactness, metallicity, internal dynamics, to¬ 
tal mass, spatial distribution and distribution of the individual 


masses, perhaps indicating differences in the Initial Mass Func¬ 
tions (IMF). It is clear that a significant improvement in the com¬ 
pleteness of the census of different associations provide key in¬ 
formation to understand star formation and evolution. 


In this context, the open cluster Messier 35 (NGC 2168, here¬ 
after M35) is very interesting in several aspects. First because it 
is a quasi-twin of the Pleiades cluster, with a similar age (Vi¬ 
dal 1973|l but a different metallicity. Second, the Kepler K2 mis¬ 


sion (Howell et al. 2014) has just provided exquisite photometric 
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time series, so that a detailed distribution of the rotational pe¬ 


riods for different masses will be available (e.g Meibom et al. 
|2009| . Moreover, M35 is not so crowded that confusion is a ma¬ 
jor problem and its size (about 1 sq. deg) is very suitable for 
follow-up multi-object spectroscopy (Barrado et al.|2001a| ). 

Several studies have provided estimates of other key clus¬ 
ter parameters such as its age, its census and its metallicity. Ta¬ 
ble [T] gives an overview of the cluster’s properties found in the 
literature in the past fifteen years. However, there is no a clear 
consensus for most of them. |Sung & Bessellj ( | 1999) 1 derived a 
distance modulus of (m - M )o = 9.60 + 0.10 corresponding to 
832 pc, a reddenin g of E(B - V) - 0,2 55 ± 0.024 mag and an 


Kalirai et al. 


(2003 pushed the cluster a lit- 


age of 200^ jqq Myr. 
tie bit further, to 912^2 pc and derived an isochrone fitting age 
of 180 Myr. More recently, [McNamara et al.| ( |2011[ ) obtained a 
distance of 732 ± 145 pc from the internal proper motions. This 
dynamical distance and the radial velocity dispersion obtained 
by |Geller et al.| ( |2010) > imply an age of 133 Myr, according to the 
same authors. An age determination based on a completely dif¬ 
ferent technique -gyrochronology- was performed by |Meibom| 
et al. ( 2009[>, yielding a value between 134 Myr and 161 Myr. 


Reimers & Koester (1988), by using the white dwarf cooling age, 
derived an age in the range 70-100 Myr, whereas Barrado ( 2001|>, 


with the same technique, estimated 180 Myr. The rotation pat¬ 
tern suggests an age slightly older than the Pleiades (175 Myr 
compared to 125 Myr). The same conclusion is reached when 
the lithium abundances are compared, since the lithium disper¬ 
sion, which is itself related to rotation for a given stellar mass, 
is no longer obvious for M35 (Barrado et al. 2001a) while it 
is very significant in the case of the Pleiades (see, for instance 
|Soderblom et al.|1993j ). 

Regarding chemical composition, [Steinhauer & Deliyannis 
( |2004| ) quoted an unpublished work with a metallicity [Fe/H] = 
-0.143 ± 0.014, very similar to the value obtained by our group 
([Fe/H] = -0.21 + 0.10) using high spectral resolution data and 
spectral synthesis, ( Barrado et al.|200~Ia| . 


The cluster is also relatively massive: Leonard & Merritt 
( |1989| ) derived a central total mass in the range 1600-3200 
Mg using dynamical models. Our own estimate, for the cen¬ 
tral of 27.5x27.5 square arcmin, is ~ 1600 M G , which allows a 
exquisite determination of a complete stellar mass function (Bar¬ 
rado et al.] |2001b so far the deepest for M35) which is totally 


compatible with the IMF of a much younger and un-evolved 
cluster such as Collinder 69 (Barrado et al. 2005] |Bayo et al.| 
201 1[), and suggests that M35 includes several tens of potential 


substellar members. 

The goal of this paper is to complete the census reaching well 
inside the substellar domain and including the cluster outskirts 
up to a radius of 1°, providing multi-epoch, multi-band photom¬ 
etry from public archives and from our own recent observations. 


2. Datasets 


We searched the CFHT, NOAO, Subaru, Kiso, SDSS (DR9, [Ahnj 
et al. |2012 ) and ING public archives for wide field images of the 
cluster. Table [2] gives an overview of the dataset found in these 
archives and in our own private archives. 

In each case, the raw data and associated calibration frames 
were downloaded. The individual raw images were processed 
using an updated version of Alambic (Vandame 2002| ), a soft¬ 
ware suite developed and optimized for the processing of large 
multi-CCD imagers, and adapted for the instruments used here. 
Alambic includes standard processing procedures such as over¬ 
scan, bias and dark subtraction for each individual readout ports 
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Fig. 1 . Estimated error on the proper motion as a function of /-band 
magnitude. The color scale represents the maximum time difference 
used for each measurement. The prediction for Gaia is over-plotted for 
comparison. 


of each CCD, flat-held correction, bad pixel masking, CCD-to- 
CCD gain harmonization in the case of multi-chip instruments, 
fringing correction (when needed for the reddest filters), de¬ 
stripping, background subtraction and non-linearity correction 
in the case of near-infrared detectors. 


Aperture, PSF and model photometry were extracted from 
the individual images using SExtractor ( |Bertin & Arnouts| 
1996) and PSFEx (Bertin [2~01 1 [ ). The individual catalogs were 
then registered and aligned on the same photometric scale using 
Scamp ( |Bertin|2006| ). 

The photometric zero-points were derived using standard 
fields obtained the same night whenever available, or all-sky cat¬ 
alogues such as 2MASS ( Skrutskie et al.|2006| l in J, H and Ks, 
SDSS in g, r, i, z and IPHAS ( |Drew et al.|2005| ) in Her. The zero- 
point rms were typically in the range between 0.02 and 0.12 mag, 
and were added quadratically to the measurement uncertainty. 
No attempt was made to calibrate the non-standard VR filters 
used with the Mosaic and SuprimeCam cameras and the He I 
narrow band filter obtained with the WFC. 


Proper motions for all the sources were computed using 
the method described in |Bouy et ah] ( |2013] ). The final residu¬ 
als on the astrometric solution show a dispersion of the order of 
11 ~ 12 mas. Figure [T] shows the estimated error on the proper 
motion as a function of /-band apparent magnitude. We note that 
the proper motion measurements are not tied to an absolute refer¬ 
ence system such as the ICRS (International Celestial Reference 
System). They are nevertheless expected to be close to the ICRS 
since extragalactic sources and background stars largely domi¬ 
nate the sample and display very little motion with respect to the 
ICRS. 

A final master catalog of all sources (338 892 in total) 
was then produced including the measured photometry in the 
grizJHKs filters. A total of 194452 sources (approximately 
half) were detected at more than two epochs and their proper 
motion are also included. Since the main goal of the survey fo¬ 
cuses on the study of M35, only objects with \/u a cos 5| <100 mas 
yr-1 and |/r,, <100 mas yr 1 were considered for the member¬ 
ship analysis. 
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Table 1 . Overview of previous studies of M35 


Reference 





Distance 

(pc) 

Age 

(Myr) 

E(B - V) 
(mag) 

[Fe/H] 

Reimers & Koester ( 

1988 

) 


70-100 



Sung & Bessell 

199 

9) 


832+39 

2oo +200 
^ -100 

0.255+0.024 


Barrado 

(200 






180 



Barrado et al 

(200k 

ii 


912+7° 

* -65 

) 

>125 


-0.21+0.10 

Kalirai et al. 

2003) 



180 



Steinhauer & 

Dehyannis 

2004 



-0.143 + 0.014 

Meibom et al 

.(2009) 



134-161 



Geller et al. ( 

2010) 



133 



McNamara et al 

(2011 


732+145 





Table 2. Instruments used in this study 


Observatory 

Instrument 

Filters 

Platescale 

[pixeR 1 ] 

Chip layout 

Chip size 

Field of view 

Epoch(s) 

Ref. 

CFHT 

UH 8K 

R, I 

07205 

4x2 

2kx4k 

29'x29' 

1996 

1 

CFHT 

CFHT 12K 

B,V,R,I,z 

07201 

6x2 

2kx4k 

42'x28' 

1999, 2000, 2001, 2002 

2 

CFHT 

MegaCam 

gT 

07187 

4x9 

2kx4k 

l°xl° 

2003, 2005,2011 

3 

CFHT 

WIRCam 

H2 

07304 

2x2 

2kx2k 

20'x20' 

2008, 2009, 2010,2011 

4 

Subaru 

Suprime-Cam 

VR-Broad 

07200 

5x2 

2kx4k 

34'x27' 

2005 

5 

INT 

WFC 

g, Ho-, He I, i, I, r, U, V 

07333 

3x1+1 

2kx4k 

34'x34' 

2000, 2002, 2003, 2005, 2008, 2009, 2010 

6 

KPNO (Mayall) 

NEWFIRM 

J, H, Ks 

07400 

2x2 

2kx2k 

28'x28' 

2008 

7 

KPNO (Mayall) 

MOSAIC 1 

VR-broad 

0726 

4x2 

2kx4k 

36'x36' 

2006, 2009, 2010,2011,2012 

8 

KISO 

2KCCD 

B,V,R,I 

1750 

lxl 

2kx2k 

50'x50' 

2002, 2003 

9 

SDSS DR9 

SDSS camera 

u, g, r, i, z 

07396 

6x5 

2kx2k 

3° 

2006 



References: (l) jMetzger et al.|11995). (2 ) jCuillandre et al.|12Q00) , (3) |Boulade et al.|12003) , (4) |Puget et al. 1(2004) , (5) |Miyazaki et al. 112002) , (6) |Ives|)1998) , (7) |Autry et al.|(2003) - (8) 
|Wolfe et al.|)2000) , (9) |Itoh et al.|)2001) , (I ()) | Ahn cl al.|)20 1 2) 


3. Membership 

We use the method described in |Sarro et ~akj ( |2Q14| > to derive 
membership probabilities of all the sources. Briefly, the method 
uses a training set - in this case the sample of candidate mem¬ 
bers from Barrado y Navascues et al.]( |20Ql| > - to identify the fea¬ 
tures that most efficiently separate the cluster’s population from 
the field and background population among all combinations of 
proper motion, luminosities and colors available. In the case of 
M35, the best set was found to include the proper motions, the 
z, H , Ks luminosities and the (r - z), (r - i) and (/ - H ) colors. 
Once the best set of colors, luminosities and proper motion is se¬ 
lected, a model of the cluster’s locus defined by the training set is 
fitted in the corresponding multi-dimensional space. A mixture 
of Gaussians is used for the proper motion diagram and princi¬ 
pal curves for the cluster’s sequence in all color-magnitude and 
color-color diagrams. Membership probabilities are then com¬ 
puted using this first model. The highest probability members 
are selected as new training set, and the models and the mem¬ 
bership probabilities are recomputed until they converge. 

This method was originally developed and optimized for 
clusters like the Pleiades, i.e clusters with a relatively large 
proper motion and photometric sequence clearly separated from 
the field and background sequences in many color-magnitude 
diagrams. Unfortunately M35 is a worst-case scenario for the 
method. Its proper motion is relatively small and its sequence 
largely overlaps with the field and background sequence in al¬ 
most all color-magnitude and color-color diagrams. To limit the 
contamination and enhance the contrast between cluster mem¬ 
bers and field and background sources, we chose to restrict the 
analysis to a radius of 1° around the cluster’s center located 
at a=06h08m54.0s, and <5=+24°20'00"(J2000). Figure [2] shows 
their spatial distribution. It shows that the north-western part of 
the survey was shallower than the rest of the survey, as fewer 
members were selected. It also shows that a number of high 
probability members have (r - i ) or (i - Ks) color apparently 


inconsistent with the cluster’s sequence. Their high probabil¬ 
ity is due to their proximity to the isochrone in all other color- 
magnitude and color-color diagrams as well as to the cluster’s 
locus in the proper motion diagram. These sources - probably 
suffering from inaccurate photometry in one of these bands - il¬ 
lustrate the robustness of the selection method based on a high 
dimensional space. 

Table ?? gives the astrometry, photometry and membership 
probability for the 338 892 unique sources detected in the cor¬ 
responding area. Note that only the sources with a proper mo¬ 
tion measurement have a membership probability (194452 in 
total), as the availability of a proper motion measurement was a 
necessary condition in our membership analysis. Using sources 
with membership probability greater or equal to 0.9, we derive 
a median proper motion for the cluster of (/ i a cos 6,fis)=( 1-99, - 
0.23)±(0.82,0.71) mas yr-1, where the uncertainties correspond 
to the scaled median absolute deviation. 

Figure |4] shows the membership distribution and cumulative 
distribution (computed from large to small) for all the sources in 
the catalog. The vast majority of the sources has a membership 
probability of «0.0. A monotonously decreasing distribution fills 
the probability range between 0 and 1.0, as the result of the lose 
constraints on the photometric and astrometric properties of the 
cluster. 

NGC 2158 is an intermediate age (1 ~2 Gyr) cluster at a dis¬ 
tance of ~3 600 pc (Carraro et ak] |2002) i located at a projected 
distance of only 25' from M35. It is therefore included in the 
area covered by the present study and possibly contaminates the 
membership analysis. To estimate the level of contamination due 
to NGC 2158 members, we select all the sources located within 
3' of its center. Given the distance and projected density pro¬ 
file of NGC 2158, these sources are expected to be in their vast 
majority members of NGC 2158. Figure [5] shows the sequence 
formed by these sources in various color-magnitude diagrams 
made of colors and luminosities used for the membership anal- 
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Fig. 3. Color magnitude and proper motion diagrams of the 194452 sources of the sample detected at more than two epochs. Left: r vs r — 
Middle: i vs i - Ks. Right: proper motion diagram. Sources with a membership probability greater or equal to 0.5 are represented in the top panels 
with a color scale proportional to their probability. The rest of the sources is represented in the bottom panels. 


Barrado+2001 
• DANCe 



93.0 92.5 92.0 91.5 

R.A. (J2000) 


Fig. 2. Spatial distribution of sources used as training set (green squares, 
from |Barrado y Navascues et al.|200 1 ( and of the 4 349 sources of the 
final sample with membership probability greater than 0.5. Background 
image: 2MASS J-band. 

ysis presented here, and compares it to the sequence formed by 
objects with a probability of membership to M35 greater than 
0.5 according to our analysis. The two sequences overlap mostly 
at the bright end, and contamination must occur mostly above 
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Fig. 4. Distribution (blue, left scale) and cumulative distribution (red, 
computed from large to small values, right scale) of membership prob¬ 
abilities for the 194452 sources of the sample 


i <17 mag. We find that only 27 sources with a probability 
of membership to M35 greater than 0.5 are located within 3' 
of NGC 2158, and 53 within 5'. Given the projected size of 
NGC 2158 core, little contamination is expected beyond 5', and 
the contamination by NGC 2158 members must not reach more 
than 1 ~2% at most. 














































H. Bouy et al.: Messier 35 (NGC 2168) DANCe 



Fig. 5. Color-magnitude diagrams of the M35 high probability members (>0.5. red dots) and of all the sources within 3' of NGC 2158 center (blue 
dots). 


4. Conclusions 

We have processed and analyzed a total of 4 867 wide-held opti¬ 
cal and near-infrared images of M35 obtained from various pub¬ 
lic astronomical archives. The images cover 18 years of time 
baseline and typically reach a depth of 21 ~22 mag in the op¬ 
tical and «19 mag in the near-infrared, while saturation arises 
around 12~13 mag at the bright end. This dataset allows us to 
derive multi-wavelength photometry from the optical to the near- 
infrared for 338 892 sources, and proper motions for 194452 
of them. We use these measurements to derive the membership 
probability of all the sources with a proper motion measurement. 
We obtain: 


- 4 349 candidate members with probability greater than 50% 

- 1 726 candidate members with probability greater than 75% 

- 305 candidate members with probability greater than 95% 


These numbers are consistent with the census presented in Bar- 
rado y Navascues et al. < 200l[ > over a smaller area. Contamina¬ 


tion is expected to be significant as the locus of the cluster in the 
proper motion, color-magnitude and color-color diagrams over¬ 
laps largely with the field and background locus. This new study 
nevertheless provides a coherent and quantitative membership 
analysis of M35 based on a large fraction of the best ground- 
based datasets obtained over the past 18 years. 
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